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Abstract 

The classical continuum mechanical model of granular media of 
rational thermodynamics results in a Coulomb- Mohr type equilibrium 
stress-strain relation. The proof is based on a two component material 
model introducing a scalar internal variable. Here we will show how 
one can get similar stress-strain relations without assuming a balance 
of substructural interactions, considering only the restrictions of the 
Second Law of thermodynamics. 

1 Introduction 

In their classical paper Goodman and Cowin derived a material model of 
porous and granular media using pure thermodynamic reasoning Ij. They 
considered a material where the density of the solid component is 7, the total 
density is p and the volume distribution function v is defined by the following 
formula 

p = z/7. 

Later the scalar internal variable u was interpreted as roughness and its (sub- 
stantial) time derivative z> as abrasion |21E1- Goodman and Cowin assumed 
a material with incompressible solid component and with a balance form 



dynamic equation for the abrasion, a balance of substructural interactions 

= V ■ h + (Tj>, 

where h is the conductive current and a^, is the production of the abrasion. 
After introducing a suitable constitutive space they investigated the require- 
ments coming from the Second Law of thermodynamics with Coleman-Noll 
procedure. Their final result was a definition of Coulomb granular material 
by the following stress function 

= (pQ - + aVu ■ Vz/ + 2az/Az/)I - 2aVz/ o Vz/ + XTr{e)l + 2/ie. 

where (3o,(3,a are material parameters, A,/i are the Lame coefficients. I is 
the second order unit tensor, o denotes the tensorial product and Tr is the 
trace. As one can see an ideal elastic behavior is coupled to the gradient 
dependent characteristic part, represented by the first two terms. 

In this paper we will show, that a similar material model can be derived 
without assuming a balance of the abrasion, with the very same weakly non- 
local extension of the configurational space as in the model of Goodman 
and Cowin. In the following Liu's theorem will play an important techni- 
cal role. Details of the different state spaces, more detailed description of 
thermodynamic concepts, the applied mathematical methods (especially Liu 
procedure) can be found in jH], regarding the weakly nonlocal extension see 
jZllHj. A thermodynamic background of continuum field theories is in jH]. 

2 Weakly nonlocal fluids - granular media 

In our treatment the basic state space of granular media is spanned by the 
density of the solid component, the volume distribution function and the 
velocity (7, i/, v). This basic state space is the simplest large deformation 
treatment and considers the possibility of changes in the topological struc- 
ture (changes of neighbours and rotation of grains) typical in fiuids. With 
this basic state we are constructing a constitutive model of a dilatant gran- 
ular material, because p is not necessarily constant. The constitutive state 
space contains gradients of the basic state variables as in case of classical 
fluids. Therefore, it is spanned by the variables (7, V7, z/, Vi^, v, Vv). The 
functions interpreted on the constitutive space are the constitutive functions. 
The space of independent variables is spanned by the next time and space 
derivatives of the constitutive variables (7, V7, V^7, z>, Vz>, V^z^, v, Vv, V^v), 
as a consequence of the entropy inequality. Here denotes the second space 
derivative. 
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In the following we will assume that the solid part of the material is 
nearly incompressible, therefore the deformation is due to the changes in the 
volume distribution function. This condition can be expressed as 

7 = (1) 

Considering ((T)) the continuity equation can be written as 

z> + z/V • V = (2) 

(m and are constraints of the independent variables, and are to be 
considered in the application of the Liu procedure. Moreover, the struc- 
ture of the constitutive state space implies that the space derivatives of the 
above equations contain terms solely from the space of independent variables, 
therefore they are also constraints. The derivative of is 

V7 = 0, (3) 

and the derivative of (j21) results in 

Vz> + Vz/V ■ V + i/VV ■ V = 0. (4) 

Finally, the balance of momentum is written as 

7Z/V + V ■ P = 0. (5) 

Here P is the pressure tensor. The requirement of nonnegativity of the 
entropy production is 

7i/s + V-j, = (T, >0, (6) 

where the specific entropy s, the conductive current of the entropy js and 
the pressure P are the constitutive quantities, functions interpreted on the 
constitutive space. With given constitutive functions the dynamic equations 
of the granular continua are ((H), Q and According to the Second Law 
we are to find these constitutive functions that the entropy production be 
nonnegative. In this way the nonnegativity will be a pure material property, 
independently of the initial conditions. Liu procedure is applied with the 
multiplier form (see [HUE]) 

ps + V ■ js - ri7 - r2V7 - r3(z> + z/V ■ v) 

- r4(Vz> + Vz/V ■ V + z^VV ■ v) - r5(7z/v + V ■ P) > 0. 

Here we introduced Lagrange-Farkas multipliers ri,r2,r3,r4 and Fs for 
the constraints ((H), ((Hj), ((21), & and ^ respectively. In deriving the Liu 
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equations one should consider that the substantial time derivative does not 
commute with the space derivative. The following identity is to be applied 

(Va) = Va - Vv • Va. 

Now the multipliers of the independent variables give the Liu equations, 
respectively. Introducing a shorthand notation for the partial derivatives of 
the constitutive quantities as e.g. = 4- we get 



pd^s 


= Ti, 


(7) 


pdy^s 


= Ts, 


(8) 


pdyS 


= Ta, 


(9) 


pd^uS 


= r4, 


(10) 


pd^s 


= prg, 


(11) 


pds/wS 


= 0, 


(12) 




- r55v.P)^ = o, 


(13) 






(14) 




- r4Z/I - Ts • 9vvP)' = 0. 


(15) 



The superscript * denotes the symmetric part of the corresponding func- 
tion. Equations (f7|)-()ll|) determine the Lagrange- Farkas multipliers. The 
solution of (fT^ results in an entropy that is independent of the gradient of 
velocity, therefore ()15|) can be integrated and a particular form of j<j is deter- 
mined. Substituting that form of into (fT!?|l - (fT^ we get two equations that 
are fulfilled if 



dy^yS = 0, (9j.v7-5 = 0, chuVuS = 0, d\/uv-/S = (16) 

The most general generalized entropy function, isotropic and second order 
in V and in V7 that satisfies the above conditions is the following 

v2 (V7)2 
s(z/, Vz/,7, V7, v) = Se{i^,'j)-m{u,-f)Y-(^i^,l) ^ ■ (17) 

Here m and a are arbitrary nonnegative functions. We can see, that 
entropy function is a concave function of the variables v and V7. Moreover, 
the entropy is independent of Vz/. Considering (fTTjl . the solution of the last 
Liu equation (fT3|) gives the entropy current as 

= -mv ■ P + ji(z/,7, v). (18) 
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Here ji is an arbitrary function. Applying (fT7|) and ()18|) the dissipation 
inequality simplifies to 



V ■ ji - V(mv) : P + p(«V7 o V7 - ud^sl) : Vv > 

If ji = and m = 1 (as usual) then the dissipation inequality can be 
transformed into a solvable form. Introducing the notation p = —vpduSe 
entropy inequality further simplifies to 

- (p + pdua^-^^)l - paV-fV-f j : Vv > 

The notation is not arbitrary, because d^Se = dpSed^p = — ^r^, where 
p is the scalar pressure according to the traditional thermostatic definition, 
corresponding to the Gibbs relation. Or, alternatively investigating a pure 
mechanical system one can recognize that our s function is still connected to 
the entropy only in some properties, therefore one can introduce the physical 
entropy s = jr. In this case p is the scalar pressure, Tq is a constant temper- 
ature and all the considerations above are valid. Let us underline here, that 
there is nothing mysterious in the above identifications that would weaken 
the above reasoning. Our whole procedure is based on the existence of a con- 
stitutive function with the nonnegative balance, a kind of general stability 
requirement. One can exploit this property and recognize the physical mean- 
ing of the corresponding quantities later. On the other hand we should not 
forget that the introduced s is a kind of generalized, nonequilibrium "physi- 
cal" entropy with rather strange variables (velocity and gradients) and there 
is no unique nomination for such quantities (coarse grained kinetic poten- 
tial?). 

The inequality above contains only the pressure tensor P as constitutive 
quantity that depends on the velocity gradient, therefore it is solvable. The 
general solution is 

P-Pe = P. = -L(Vv), 

where P„ is the viscous part of the pressure, L is a nonnegative constitutive 
function. Here we assumed a symmetric pressure tensor as it is usual for 
materials without internal moment of momentum. Furthermore, the equilib- 
rium, static part of the pressure is Pg and it is defined as 

Pe = (p + pd,a^^)l - paV-fV-f (19) 

A material defined by this equilibrium pressure tensor will be called 
Coulomb-Mohr material. 
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3 Coulomb- Mohr materials 



The pressure tensor fll9|l has several remarkable properties. First of all, the 
static shear stress is not zero, and V7 is an eigenvector of Pg. Furthermore, 
introducing an arbitrary direction with a unit vector n, one can define a 
pressure normal to the direction n as 

N :=n Pf. n = p — pa(V'y ■ n)^, 

where p = p + pdyaiV'^Y' 12. Denoting the shear pressure hy S := P — N 
one can get 

S'^ + N^ = (P, . n) ■ (Pe ■ n) = p2 _ 4pap(V7 " + Kp<^f{^lf{^l " 
After a short calculations follows, that 

S"^ + (A^ — tY = where t = p ~ r and r = pa(V7)^. 

Hence, the possible Mohr circles in the material have a special structure, 
their envelope from above is a straight line, the material satisfies a failure 
criteria of Coulomb-Mohr kind, as the Coulomb material of Goodman and 
Cowin. One can check that this is really a failure criteria, the second deriva- 
tive of the entropy function (fTTj) become semidefinite on the Coulomb-Mohr 
line and the line separates the regions of the state space where the thermo- 
dynamic stability is violated from the region where it is fulfilled. 

Material stability in mechanics is far more complex matter than in equi- 
librium thermodynamics of gases and fiuids where it is connected to the 
appearance of new phases E] • On the other hand in a thermodynamic 
approach the loss of thermodynamic stability means the violation of dy- 
namic stability (e.g. the corresponding partial differential equations loose 
their hyperbolicity). Stability losses of this kind where treated in connection 
of rigid, microcracked materials in [13j. There the microstructure was char- 
acterized by a vectorial dynamic variable D (representing the average crack 
length, a damage) and two general free energies of the microcracked materi- 
als were introduced. There with special material parameters one can arrive 
to a Coulomb-Mohr failure criteria and pressure tensors similar to (fTI?|) . but 
with a local variable called damage vector D, instead of the gradient variable 
V7. The role of V7 is similar to D in every sense. Recalling the quadratic 
form of the entropy function the loss of thermodynamic stability is not a di- 
rect consequence of an increasing IV7I. The loss of stability is more involved 
and always connected to changes in the stress/deformation state and appears 
only in perpendicular to V7, in a 'shear' direction. The complexity of the 
different kind stability losses is investigated e.g. in [Ti| I15j. 
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Figure 1: The envelope of the Coulomb- Mohr circles as a boundary of the 
stability 

Our basic state space and configuration space was that of a fluid. In the 
one component case, when the configuration state space is (p, Vp, v, Vv) an 
ideal Euler fluid and the viscous Navier-Stokes fluids appear after similar 
derivations. Remarkable, that a higher order weakly nonlocal one compo- 
nent fluid, where the constitutive state space contains also the second order 
derivative of the density, results in the Madelung fluid (known from the hy- 
drodynamic version of quantum mechanics) in the ideal, non dissipative case 



According to the comparison of Kirchner and Hutter the model of Good- 
man and Cowin is the most robust one among several local continuum gran- 
ular material models of porous and granular media (with scalar internal vari- 
ables). Interestingly, the other compared models are all can be considered 
as relocalized (according to the terminology of jZj), because they can in- 
troduce nonlocality through a generalized entropy current. A different con- 
tinuum granular model of Aranson and Tsimring introduces also a scalar 
internal variable with a Ginzburg-Landau dynamics (assumed in an ad-hoc 
way) ^3 QHl Onj- Their starting point is the granular kinetic theory and 
thorough calculations and simulations try to understand the connections in 
that direction. Therefore they model does not compatible to a realistic static 
case. In understanding the granular phenomena and especially to create an 
applicable continuum model is a great challenge of physics where definitely 
a large extension of our understanding of continuum concepts and therefore 
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the extension of continuum theories is necessary |20j. The simple model 
shown in this work is hopefully is a step in that direction demonstrating the 
capabilities and simplicity of thorough constitutive reasoning. 

Finally, it is important to observe the differences in the conditions of the 
derivation and also in the final results of the model of Goodman-Cowin and 
the present one. Here the fluid like state space and the large differences in 
the compressibility were essential conditions in the derivation. There was 
no need to introduce a balance of substructural interactions (but of course 
the state space can be extended to include 7 and one can investigate the 
consequences). We have got a Coulomb- Mohr failure criteria in the static 
limit in both cases. However,the structure of the equilibrium pressure tensors 
and therefore the implied physical reasons were completely different. To 
investigate the physical relevance of the two models further work is necessary. 
In this direction are very important the known solutions and generalization 
of the Goodman-Cowin model like the flow calculations of fH] I22j. Another 
possibility to test the differences could be to look into static failure data of 
Coulomb-Mohr materials. Related experimental data regarding porous rocks 
is published recently by Vasarhelyi [231 1211 12S] 
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